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The physics of the Wilson line leads to new developments in high temperature 
particle physics. The main tool is the effective action for a given fixed value of 
the phase of the Wilson line. It furnishes a gauge invariant infrared cut off, and 
yields for small values of the phases a systematic procedure for obtaining a power 
series in the coupling g and glog(l/g). It breaks the centergroup symmetry of the 
gauge group only at high temperature so leads to domain walls disappearing at low 
temperatures. It shows long lived metastable states in the standard model, SU(5), 
SO(10) and its SUSY partners, with possibilities for CP violation and thermal 
inflation. 



1 Introduction 

Phase transitions have been a recurrent theme at this meeting. We have heard 
about the quantitative progress in understanding the electro-weak transitionEl 
and along similar lines the SU(5) GUT phase transitiona and to some extent 
the deconfinement transitionEl. 

Despite their differences there is a common feature: all of them have an 
order parameter, the Wilson line. The Wilson line is strictly speaking an order 
parameter, when there are no complex representations of the gauge group 
present in the particle content. Only then it will be strictly zero in one phase 
and non-zero in the other. 

Usually the Wilson line is mentioned in discussions of the quark gluon tran- 
sition (apart from the chiral condensate), not of the other transitions, like the 
electro-weak and GUT transitions where the Higgs fields are the protagonists 
in the transition. 
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In this paper we will discuss the role of the loop; in particular its phase, 
and the dependence of the free energy on that phase. The latter is called the 
effective action and is the main mathematical tool. 

For small values of the phases this effective action is a gauge invariant infra 
red cut off version of the free energy and the Debije mass. As such it is a nice 
starting point for the perturbative evaluation of these quantities □. Infra red 
singularities upset the usual series in g 2 , and their computation necessitates a 
cutoff. This is briefly discussed in section |[ 

For values of the phase in the centergroup one finds metastable minima. 
They have one common property: they are very long lived on any cosmological 
scale. This is true for the standard model and for GUT models like SU(5) 
and SO(10) and their SUSY versions. A rather detailed description of where 
the minima are and how they decay is found in section ||. Then we turn to 
the other order parameters, the Higgs fields. The potential has qualitatively 
different behaviour due to the Wilson line condensate, in case the Higgs field 
carries Z(N) charge. This is the subject of section |[ 

In section 5.3 we turn to the physics of these metastable states. As the 
CP violating properties of some of the metastable states have been analysed 
elsewhereH, we limit ourselves to thermal inflation. This happens in the SO(10) 
model and its SUSY analogue. 

The question of how the universe can arrive, in such a metastable state 
arises naturally B. There has been debate on thisB , as well as on the thermo- 
dynamical properties, of the metastable states and the physical relevance of 
the phase of the loopB. On the latter two the reader can find partial satisfaction 
in section 0. 



2 Wilson line as order parameter 

In this section we will resume some well known facts about gauge systems at 
high temperature in equilibrium. First the academic but instructive case of 
pure SU(N) theory is reviewed. Then we introduce the Wilson line as the 
order parameter that describes walls, and the effective action that controlls 
its behaviour. In the last subsection we discuss the boundary conditions that 
trigger localised walls. 



2.1 Pure gauge theory 

Let us consider a pure SU(N) gauge theory i. 

Take an elongated box in 3D with size L 2 r L z . The size in the z-direction is 
by far the largest in order to have walls that separate vacuum states. Boundary 
conditions on the gauge potentials are periodic. 
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Physical states are by definition invariant under periodic gauge transfor- 
mations. This means Gauss' law is satisfied everywhere. Consider now a gauge 
transformation periodic in the transverse directions, but periodic modulo a 
centergroup element expi^M m the z-direction. Such transformations Uk dif- 
fer from one another by a periodic transformation, so have all the same effect 
on a physical state. Like a periodic transformation they do commute with 
the Hamiltonian. So we can diagonalise them simultaneously. Acting on the 
physical states they have eigenvalues all of which are again Z(N) phases|| 

We can easily get a physical picture of what these states are: start from a 
state on which all operators Uk have eigenvalue one. This state has no electric 
flux, only glueballs are thermally excited. Now take a string running from z = 
to z = L z given by the path-ordered operator TrPexpz J dzA z (x,y, z) acting 
on the no flux state. This state has one electric flux, since the eigenvalue of 
Uk is expik^-. If we add another string k will be replaced by 2k in the phase. 
Note that a state with a given flux e can have e mod N strings. 

We can create strings independently in all directions. This means that 
a state will by given by a flux vector e. Such a flux state is obtained by 
acting with a projector Pg on a physical state. This projector is related to the 
operators Us by a Z(N) Fourier transform: 

1 - 2tt 

k 

Since flux is conserved (the Hamiltonian commutes with the Us), we can 
define a flux free energy Fg by the Gibbs trace over flux e states: 

ex P^-r = ^<e| exp -^-|e) (2) 

These flux free energies should tell us how the strings behave. Now we expect 
at low T confinement, hence a string tension p. Let us take the case of one 
flux in the z-direction. Then at low T (<C p): 

Fi - F = P L Z (3) 

For high enough T(T C ) this behaviour changes into: 

F 1 -F = L 2 exp^P (4) 

The physics of these equations is simple. At low T creation of one string has 
a small probability ~ exp— ^L z . At the same time the difference between 

"That's because is periodic, hence does not change the physical state. 
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the free energies increases as the length of the box increases. The unique 
groundstate is given by Fg. 

As the temperature goes up so goes the probability for exciting strings. 
The number of strings present in the Gibbs sum grows and so does the en- 
tropy. At T c the entropy overtakes the energy p and above T c the free energies 
become exponentially degenerate. This means that the Z(N) symmetry is 
spontaneously brokcn.j] 

The question is now: does the parameter a that controls the decay of the 
flux free energy correspond to a surface tension between two regions degener- 
ate in energy? If so, there must be an order parameter telling the difference 
between the two degenerate states. 

2.2 Domainwalls as Wilson line profiles 

To study domain walls one has to introduce the duals of the flux free 
energies. They are called "twisted" transition elements and are defined by 

z k = Tr phys exp {-H/T)U % (5) 

The Z^ are related to the Fg by the formula |l| for the electric projectors, 
that is, by a Z(N) Fourier transform. 

They can be rewritten into 4D "twisted" path integrals. To see how this 
works, take Ut = 1. This corresponds to the well-known pathintegral with an 
integration over periodic Aq coming from the physicality constraints on the 
states. 

The presence of Uj: does change the state the path integral starts with at 
time r = 0: it creates a center group discontinuity when going from one side 
(z = 0) to the other (z = L z ). Then we go in the time direction to the point 
t = 1 /T, z = L z , keeping this discontinuity. Going first in the time- and then 
in the z direction one meets no discontinuity, so we have a vortex in the z — r 
planes with strength expifc^. Consider a path ordered Wilson line 

P(A (x)) = J P{exp^ A a (r,x)dr} (6) 
Jo 

''This behaviour has been foustdi analytically in gauge Potts models and in gauge the- 
ory by Monte Carla-simulations til- For other mechanisms of symmetry breaking at high 
temperature see ref£3 

c Strictly speaking, Aq need not be periodic in the Euclidean time direction. Constraining 
it to be periodic does not affect the thermodynamical properties 
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at z = 0. When we push the line to z — L z it will pick up the phase of the 
vortex, as soon as it crosses the center of the vortex^ 

This is an important property of the line: at high T it will have a non-zero 
expectation value and change its phase going from one side of the box to the 
other. This answers the question we posed at the end of the last paragraph: 
the parameter a corresponds to a surface tension, the surface separating two 
phases where the Wilson line takes different phases. Note that we have not as 
yet a three dimensional interpretation of the Wilson line. As it stands it is a 
Euclidean path integral object. 

The twisted transition element has been analysed numerically in its four 
dimensional path integral form. EH In three dimensions also the profiles p(z) 
of the Wilsonlines have been measured, clearly indicating domain walls above 



2.3 Wilsonline and heavy quarksource 

We would like to associate the Wilsonline in the periodic system with the pres- 
ence of a heavy quark. This in order to have a three dimensional interpretation 
of the line and of the twisted transition element. 

However in a periodic volume with Gauss' law everywhere imposed a single 
heavy quark source cannot exist. So we have to drop the periodicity and search 
for convenient boundary conditions, that replace our twisted box and introduce 
a localised wall. In this article we will not explain this in deiail, but the strategy 
can be read off from a simple Z(2) lattice gauge model. EJ The model is the 
limit of an SU(2) gauge model where all particles are made very massive by 
adjoint Higgs multiplets. This leaves only Z(2) centergroup transformations, 
as they still commute with the adjoint Higgses. There one can show that at 
high temperature a domain wall is created with a profile of positive energy 
density. The rigour is the same as used in proving that the Ising model can 
produce walls. 

3 Effective action, perturbative expansion, and infrared cut off 

The surface tension a was computed in perturbation theory. Also the profile 
of the Wilson line. The computation leads in a natural way±o the realisation 
that the loop serves as a gauge invariant infra red cut-off □ in the effective 
action. This is the subject of this section. 

d Clearly such a line in the z-direction will have this property too. However, it will not 
havp-p VEV at high T, contrary to the Wilson line. For a more complete discussion, see 




ref. 
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The effective action U for the Wilson line is defined for any given profile 
p(z) as follows: 



exp-^M^= J DAiiS (p- P(A ))exp-±S(A) (7) 

To avoid clutter in eq. we left out the z dependence. P(Aq) stands for 
the normalised average over the transverse directions of the Wilson line. Hence 
it need not be unitary. 

The effective action gives the twisted transition element (||) by integrating 
over all profile configurations with boundary conditions appropriate to the 
twist. In the time direction the statistics of the boson fields imposes periodicity, 
for fermions anti periodicity. 

Let us note that the Wilson line @ transforms under periodic gauge trans- 
formations as an adjoint. So only the phases are gauge invariant. Hence the 
loop P in eq. stands for TrP, and to get all eigenvalues one should admit 
winding the loop several times (from 1 to N-l) before taking the trace. 

U(p) has a symmetry under gaugetransformations that are periodic in the 
time direction modulo a centergroup clement. Such a transformation indeed 
leaves the action and the measure invariant, but changes the Wilson line by 
that same phase. Hence we have: 

U(pexpik^r) = U(p) (8) 

The same stays true when we admit fields with no Z(N) charge. But fields that 
do carry Z(N) charge will see their boundary conditions in the time direction 
changed by the phase. Hence the relation (|^) is no longer true; it only holds 
approximately, the better the larger the masses of those fields are. 

In what follows we consider the perturbative expansion of U. This implies 
we work at very high T, where the VEV of the Wilsonline has modulus one, 
so only phases: 

p(z) = expi27rC(z) (9) 

The matrix C is in the Cartan subspace of the Lie algebra of SU(N) and 
contains all the eigenvalues Ci,i — 1....N, with the constraint that they add 
up to zero. So our perturbative expansion will be around C as background. 

U(C), the effective action, will take the general form of a kinetic and a 
potential part: 

U(C) = I dz \K{C)Tt{2tiT^-C) 2 +ir 2 T 4 V(C)j (10) 
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Both K and V are Z(N) periodic just as U(C) in eq. [|, and will loose that 
property when there are Z(N) charged fields. 

At large values of L tr the pathintegral over all profiles is dominated by the 
extremum of U{p). At large values of \z\ the wall will be determined by the 
behaviour of kinetic and potential terms at C=0. The behaviour of the wall in 
that region is determined by the large distance behaviour of the correlation of 
the phase of the loop. Hence it will be dominated by the Debije mass, not by 
the magnetic glueball mass, since the correlation is odd under CTLL3. 

3.1 Background field expansion 

The expansion introduces fluctuation fields Q M around the background C: 

A ll = 2TTTC8 ll , + gQ lx (11) 
Now we expand the Wilson line around C . This gives: 

P(A ) = e*pi2irC + I dTexp{i2irTCT)igQ a (T)exp(i2TrTC(l/T-T)) 

i i 
r t r ~ 

+ / dT\ I g?T2 exp (i2TrTCTi)igQo(Ti) exp (i2i:TC(t2 — ti)) 

JO Jti 

igQ (r 2 ) exp (i2nTC{l/T - r 2 )) + 0{g 3 ) (12) 

This expansion satisfies the constraint in the effective action up to terms 
of order g. The delta function is then expanded around these terms. Since C 
is diagonal the trace of the order g term tells us that static diagonal Qo should 
not be integrated over. On the other hand expanding the action gives a linear 
term containing only static diagonal Qq. So C is the correct minimum for the 
expansion. _ . 

We have to choose gauge fixing and will take it of the form til: 

S Bf = -Tr{^D a Q Q + d.Qif (13) 

which reduces to usual background field gauge fixing for £' = 1. £' — > gives 
us static background gauge, in which Q decouples from the loop in eq. \L2. 

So in general thermal fluctuations of the loim expectation value will be 
0(g 2 ), and in static background gauge absent. EZl This is important for the 
tunneling results of next section. The fluctuations will contribute to the effec- 
tive action, and render it gauge independent. E3 

Another remark concerns the propagators of the quantum fields Q: with 
our gauge choice the background field will enter the quadratic part of the action 
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through the covariant time derivative: TrQ ^.(—Dq — 3 2 ).Q Al . The background 
enters only through commutators. Let us denote the field with all components 
zero, except the one on the row i and column j by Q tJ . Then the inverse 
propagator will look for all polarisations like: 

({2irT) 2 (n + C z -C.f+pf (14) 

This shows that the static configurations (n = 0) are screened by the 
phases, except where the differences are integer valued. This happens precisely 
in centergroup values expi2nC — expifc^. Fields diagonal in colour are not 
screened. From those diagonal fields the constraint in eq. |7] eliminates part. In 
other words the [/(l)^ -1 subgroup (including for the SM QED) keeps a infra 
red problem. 

The same is true for any particle species that has no Z(N) charge, e.g. is in 
the adjoint representation. Actually for those species all centergroup elements 
look the same. This confirms on the perturbative level what we found quite 
generally in the previous subsection. 

For particles in the fundamental N-dimensional representation the inverse 
propagators look like: 

((2 7 rT) 2 (n + C J ) 2 +p) 2 (15) 

for the i-th component. Now there is still screening in the centergroup 
elements. 



3.2 Walls and profiles in perturbation theory 

The effective action is the key quantity to compute. In the limit of very large 
L tr its extremum is going to dominate the path integral over profiles. For the 
one dimensional action in eq. ^ this extremum is simply given by the equations 
of motion for the profile: 

d 



JK(C)2itT—C - ttT 2 JV(C) = (16) 
oz 

The eigenvalues G\ = —Ci in SU(2) are parametrised by 2C = C\, V(C) is 
periodic mod 1 in C, as is obvious from the discussion on Z(N) symmetry in 
the previous sections. It is normalised to zero for C=0. Solving this equation 
with the boundary conditions C=0 and C=l (p=l and p=-l) at z — ±oo gives 
a profile, that is controlled at large \z\ according to the general arguments just 



above subsection |3T| by the Debije mass 

m D = 77 hrn * / \„i (17) 
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Let us see how this works in perturbation theory. 

One can. readily do this in the background expansion and one finds to two 
loop orderlij: 



K(C) = l/g 2 (l+g 2 K 2 (C)) 
V(C) = C 2 (l - \C\) 2 - 5|^C 2 (1 - \C\) 2 (18) 

The kinetic term K%{G) behaves well for all C, except at small values, 
where it develops a pole. This is in contrast with our general expectation of 
the large z behaviour. It should be governed by the Debije mass in eq. |l7|. 
This equation, together with eq. [l^, tells us that the correction to the lowest 
order result ^N/igT for mo is (1 + 0{^))\ However the work of RebhanS 
shows that this correction is not 0(g 2 ). It is larger, of order <?log(l/<?) ! 



3.3 Postmortem 

Our seemingly disastrous result, the pole in the kinetic term, just indicates, 
that for small C, naive perturbation theory breaks down. For small values of 
C, say 0(g 2 ), the screening of the propagators in eq. |lj becomes ineffective. 
For those values of the background one has to do first the integration over all 
non-static configurations, and obtain an effective 3D action 53. citeginsparg 
This integration will induce terms of order g 2 (the one containing the lowest 
order Debye mass term for the static fourth component of the quantum field) 
and higher. 

The new action will only depend on the static quantum fields, the back- 
ground C and inherits the gauge fixing in eq. |l3|. It will receive contributions 
from the non-static part of the constraint. The latter is gauge dependent and 
absent in static background gauge. Like in the 4D calculation they are the 
guarantee for gauge choice independence. n 

We have calculated the Debije mass from this effective action. □ In Feyn- 
man background gauge there is only one 3D diagram contributing. The result 
coincides with that of Rebhan. 



4 Effective potential for the Standard Model and beyond 

In this section the one and two loop potentials for the Standard model and 
beyond are analysed. 

There are two interesting types of minima: the absolute minimum and the 
metastable minima. Ell It will turn out that the latter are very long lived on 
cosmological scales. 
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Let us first analyse the contribution to the potential of a generic particle 
species coupled through the covariant time derivative to its Wilson line phase 

X s • 

D (x) = 2irT(n + x) (19) 
The one loop result for a complex boson neglecting its mass is then: 

V h {x) = Tr T D 2 (x)-Tr T=0 D 2 (x) 

= i^ + ^Wd-Nf (20) 

The trace at temperature T stands for the sum over Matsubara frequencies 
and a regulated integral over 3 — 2e dimensions. We subtracted the T=0 
contribution, and get, because of eq. [l^ a result which is periodic in x mod 1. 
The contribution of a Dirac helicity doublet is simply: 

Vf {x) = -V b {x+\) (21) 

The minus sign comes from the fermion determinant and the shift from 
the anti-periodicity for fermions used in eq. |l9[ 

Note that the sum of the two is antisymmetric when x is shifted over 
\. So SUSY theories, though broken at high temperature, have this discrete 
symmetry for every species. 

Now we have to specify what x is in terms of the various Cartan group 
charges of our gauge group. 

For the SM and SU(5) we are working in a four dimensional space of 
phases. So our matrix C (eq. [|) can be conveniently described by an SU(5) 
matrix in which the phases for colour SU(3) are given by q and r, of weak 
SUl(2) by t, and of weak hypercharge U(l) by s: 

,qrsqrs2qststs-. , „. 

diaq(- + , - , ,- + -, h-) 22 

yv 3 2 3 3 2 3 3 3 2 2' 2 2> v ' 

So e.g. for the righthanded electron x equals —\s. 

All what is left to do is to sum over all particle species. We then normalise 
to zero by subtracting out the Planck free energies in eq. |2(] and eq. ||l]. The 
result is plotted for SU(3) in fig.l, and for the SM and SU(5) against the weak 
hypercharge phase in fig. 2. The Higgs content of SU(5) was taken to be the 
5 and the 24. For SO(10) the 16 and and 45. For SUSY SO(10) we took the 
complexified 45 and 54 to guarantee a good SU(5) limit. 

For SO(10) the plot in fig. 2 is in terms of a U(l) charge u, orthogonal to 
the SU(5) charges. The reason is that SO(10) contains SU(5)xU(l). 
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Figure 1: 5(7(3) potentials on the elementary cell (without and with six massless fermions). 



All minima are in the centergroup elements. The trivial centergroup el- 
ement contains the absolute minimum. The other centergroup elements are 
metastable points ( though not all, as in the SM). There would have been 
degeneracy, had all the species been Z(N) neutral. But the fermions, the fun- 
damental Higgs have Z(N) charge, and lift the degeneracy. 

How can we be sure of seeing all the relevant minima? This will be ex- 



plained in subsection 4.1 



Then there are the tunneling rates of th ese m etastable states into the stable 
state. They will be discussed in subsection 4.2. Lastly the thermodynamics is 
discussed. 



4-1 Centergroup lattice of SU(N) and SO(2N) groups 

We begin with the lattice of points in the Cartan subspace, that contains all 
centergroup elements. Such a lattice is easy to generalise from the the case of 
SU(3) in fig.l. 

Consider the set of N basis vectors ex, e%, ejv given by the diagonal NxN 

matrices e\ = ^diag(l, 1, , 1,1 — N), e-i — j^diag(l — N, 1, ....1), = 

j?diag(l, 1, 1 — N, 1). Their sum adds up to zero. Taking any linear com- 
bination L of the ek with integer coefficients will give a centergroup element 
exp ilixL. The inverse is true too: all centergroup elements are to be found on 
this lattice. 

This lattice contains many elementary cells on which the potential is iden- 
tical. This is because of the symmetries in eq. |l9|. 

A convenient cell is starting in 0, and formed by the succession ei,ei + 
e%, , ei + &2 + + &N-X- Inside this cell there are conjugated points ei 
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Figure 2: Potentials in the weak hypcrchargc direction s (u for 50(10)). 



and ei +e2 + .... + ejv— 1, ei + e2 and ei +e2 + .... + cat_ 2 , etc. related by charge 
(or CP) conjugation, and with complex conjugate values for the centergroup 
elements. 

The cell of SO(10) is related to that of SU(5) because SO(10) or rather its 
covering group Spin(lO) contains SU(5)xU(l). The centergroup is Z(4) as one 
can easily check from the 16 dimensional spin representation. We normalise 
the u variable by fixing P=-l at u=l and all other phases zero. Consider the 
cell of SU(5). Shift from ei by — ^ in the u direction to get the vector u\ with 
P=-l. Shift ei + e 2 by in the u direction to get u 2 with P=i. The complex 
conjugate of u 2 is u 3 and is obtained from ei + e 2 + e 3 by shifting over — 
and that of u\ from — e 5 by shifting over i. This generalizes easily to SO(2N). 
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Where are the minima of the potential? If only Z(N) neutral fields are 
around there is a proof to all orders in perturbation theory that the minima 
are in the centergroup elements E3. 

As discussed above, fermions and Higgs fields lift the degeneracy. We have 
verified numerically that no other metastable points develop. So we have the 
following working hypothesis: 

Any metastable or stable minimum must be in the centergroup. 

A last question concerns how the minima of say SU(5) do give rise to the 
minima in the standard model, when decoupling the heavies in SU(5). 

Our working hypothesis will prove useful here. The heavies in SU(5) are 
given a common mass m, that we switch on from its value in SU(5) to oo 
in the SM. They are all multiplets of the SM gauge group, so switching on m 
will cause a movement of the minima of SU(5) along the centergroup elements 
of the SM, to wit Z(3)xZ(2)xU(l). Since the process is continuous the minima 
can only slide along the six U(l) lines corresponding to the discrete Z(3)xZ(2) 
group. This is shown in fig. 3 for the U(l) line q = r = t = 0. 

SU (5) 




Free Energy 



Figure 3: Potential interpolating between SU(5) and SM. 
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Similarly the evolution of the SO(10) minima into SU(5) minima occurs 
along the 5 U(l) lines, that we used in the construction of the SO(10) cell out 
of the SU(5) cell. 



4-2 Tunneling paths and tunneling rates 

Semi classical methods are appropriate, in our situation with small coupling, 
to compute the tunneling rate T per unit volume: 

r = aT 4 exp-S b (23) 

The dimensionless quantity a derives from the determinant of the fluctu- 
ations around the. classical bounce configuration B (for QCD and the SM this 
was done in ref.Ea) 

The computation of the bounce consists of finding the solutions to the 
Euclidean equations of motion. The temperature T is much larger than the 
scale of the 4D Euclidean bounce, we can do the calculation in 3DcJ for a 
radially symmetric bounce: 

d 2 B 2dB ,, s 

+ - =-V'(B) (24) 
dr r dr 

Boundary conditions are lim r ^oo = B ms and 4^r — at the starting point fn. 
Our bounce is a multicomponent object (4 or 5) and B ms is the metastable 
point in the elementary cell. We study the decay to the stable vacuum at the 
origin, so ro is near zero, depending on the thickness of the wall. We have 
developed a numerical method that solves the equations of motion and can 
handle the multicomponent field. 

Luckily there is a lot of symmetry in our problem, that allows an educated 
guess of the actual tunneling path. So we can reduce the problem to a one 
component bounce, solvable with the undershoot-overshoot method. 

The bounce path from a given metastable point in the cell is simply given 
by the straight path from that point. The direction of the bounce path is 
always inside a valley of the potential except for the SM where it stay near to 
a valley. This was corroborated by our numerical method for the simple groups. 
For the SM the deviation from a straight path was small but significant ||. The 
results for the bounce actions are presented in table 1. 

What is striking is the large value of all bounce actions. Tunneling tem- 
peratures T tU n can be estimated by the following simple- minded argument. 
Tunneling occurs at a time t c determined by Tt 4 ~ 1. Matter in our metastable 



Not all metastable minima in fig. 2 lie in the same cell as one can easily check. 
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Table 1: Numerical data on the bounces. 





metastable points 


Bounce action 


Critical radius 

(R/gT) 


SM 


s= 1.28, q = r = t = 


6.99 


2.7 


q = l,a = -0.25, r= t = 


2.64 


2.0 


MSSM 


s = 1.761,3= r = t = 


37.60 


3.25 


SU(5) 


p = exp (i47r/5) 


10.63 


2.3 


p = exp («27r/5) 


86.63 


5.3 


susy SU(5) 


p = exp (i47r/5) 


9.92 


1.75 


p = exp (i27r/5) 


50.47 


3.3 


SO(10) 


p = -l 


20.99 


2.75 


p = i 


308.8 


7.4 


susy SO(10) 


p=-l 


36.87 


2.2 


p = i 


286.07 


5.1 



states has qualitatively the same free energy as in the stable vacuum. Only the 
pressure — ir 2 T 4 V(C ms ) is lower in the metastable states. So we will assume 
the relation t ~ M Planck /T 2 . Then T tun — M Planck exp -S b /A , using eq. ||to 
eliminate t c . Putting in the numbers from the table shows that the metastable 
states easily survive all known cosmological transitions. 

It means that the behaviour of the potential as function of the Higgs fields 
will be very relevant for how the metastable states will decay. This is the 
subject of section ||. 

Fig. 4 shows the effect of the two loop contribution, the couplings were 
taken to be a — (3 = Ai = A2 = 0.1 and g = 0.4. The sign of the contribution-Is 
uniformly negative. That this was the case at C — was known since longEi 



4-3 Thermodynamics 

The calculation of the surface tension El motivated a renewed interest in the 
thermodynamics of the minima. In all the models considered in the previous 
subsections the energy and entropy densities of the minima are positive. But 
the thermal boundary conditions on fields with non-trivial Z(N) charge are 
changed in the metastable minima by a Z(N) phase, and as a consequence the 



15 



Figure 4: 5(7(5) potential: one-loop (dashed) and two-loop (solid). 



occupation number is changed i: 



= ,E I ■,2^^ 1 ( 25 ) 

exp (f + ikjf) ± 1 

So in the case of fcrmions one finds a purely imaginary value for the fcrmion 
number. So far no reasonable interpretation has been found for this behaviour. 
States with complex conjugate Z(N) phase are related by CP conjugation. Our 
phases have the C and P transformation properties of a diagonal Aq, so are 
CP odd. This singles out the self conjugate states like the one in SO(10) with 
P = — 1. There, fermion number is zero. 

Let us finally mention a case without any thermodynamic anomalies: large 
N pure gauge theory. To see this, consider the free energy without the Planck 
contribution: 

V(C) = ^T i Yl C?.(l-|Cy|) 2 (26) 

where = C{ — Cj . Due to the permutation symmetry of the variables Ci , 
this potential has an extremum in the barycenter of the elementary cell, in 
the one loop approximation this is an absolute maximum. Let us compute its 
value V max . In the barycenter the value of Cy is Kf- ■ Substituting in eq. 



26 
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gives: 

1 



V max = —(N< 1) ^1 + O(^) j (27) 

To leading order we find precisely the Planck free energy, but with the 
opposite sign! Hence adding the Planck free energy we have positive energy 
and entropy density, except in the maximum where both are zero. Hris result 
remains true in the next order, since it gives positive contribution. 

For finite N there is a small cap around the maximum, where both are 
negative. 



5 Combining Higgs and Wilson line potential 

In this section we will compute the classical and part of the one loop con- 
tribut ion to the mass term in the Higgs- Wilson line potential. Already on the 
classical level a few subtleties merit attention. 



5.1 Classical effective action 

It is best to illustrate the problem by the specific examples of the adjoint Higgs 
field £ and the fundamental Higgs field H . 

The terms of interest in the original 4D action are the kinetic terms of the 
two Higgses, minimally coupled to the gauge fields: 

S = J dxdr[Tr(D^) 2 + (D^HyD^H + ....] (28) 

To get the 3D effective action one would guess that the classical contribu- 
tion would just be the 3D reduction of the kinetic terms. That would give us 
terms like Tr[C, S] 2 . 

This cannoi be correct. We know already on general grounds that the 
effective action for Z(N) neutral fields like £ is identical in all centergroup 
elements. But this commutator term is not. 

What couples in the 3D action to the Higgs is the phase of the loop p(C). 
We define the matrix logp(C) through the diagonal form of the argument, and 
under a gauge transform Q we have 

\ogp(C n ) = mogp(C)& (29) 

Then the coupling of the Higgses to the Wilson line is given in a gauge 
invariant form in the effective action U{C) by: 

U(C) = J dzT 2 [Tr[Im\ogp(C),J:} 2 + (Im\ogp(C)Hy(Im\ogp(C)H) + .... 

(30) 
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The VEV of the Wilson line induces mass terms that merit comment: 

i) Thc adjoint mass term disappears in the centergroup, since the commu- 
tator vanishes there. 

ii) The mass term for the fundamental Higgs does not vanish in the center- 
group. Define |91ogp(C)| < n, then the mass term is the same for conjugate 
vacua. 

iii) Outside the centergroup the induced Higgs masses are not necessarily 
the same, i.e. SU(5) invariant. In general, the potential is not anymore function 
of only TrY? etc., but contains other invariants involving p(C). 

5.2 Quantumcorrections to the Higgs potential 

The quantum corrections are crucial to the understanding of the temperature 
dependence of the Higgs potential and hence for the occurrence of a phase 
transition. These corrections have been calculated for the trivial Wilson line 
condensate. 

In the presence of a non-trivial condensate these corrections will be de- 
scribed below. 

Let us take the concrete example of SU(5). The Higgs potential in the 
condensate characterised by expifc^p takes the form: 

tt 2 T 4 V = M 2 (r,fc)TrS 2 + m 2 (r,fc)i7 t i/ + Ai(rrS 2 ) 2 + A 2 rrS 4 

+aH^HTrY? + f3H^ 2 H (31) 

The mass terms are temperature and k dependent through loop correc- 
tions. In the one loop calculation that we completed the £ mass has the 
following corrections in Landau gauge: 

(188 \ T 2 

104Ai + — A 2 + 60.g 2 + 20a + 40 J- (32) 

M 2 (T, 1) = (W + ^A 2 + 60.g 2 - fa - §/}) ^ (33) 

The mass term contains a dependence on the condensate uniquely through 
the couplings to the fundamental Higgs. This is a consequence of the £ prop- 
agators not depending on the condensate. The induced terms will not change 
much the Higgs potential, so the SU(5) transition in a non-trivial condensate 
will not change essential features. 
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On the other hand the transition in SO(10) and the electro weak are 
changed by induced mass terms of order one: 



m 2 (T) = -m A + T 2 (Imlogp(C)y + O(couplings) (34) 

5.3 Some physical consequences 

Striking features of the potential are: 

i) CP conjugated states in the SM. 

ii) Occurrence of thermal mass terms of order unity for Higgs fields with 
non-trivial Z(N) charge. _ 

The first point has been discussedtl and we will not go into it. they involve 
pe se states with imaginary fermion number. 

The second point has two important consequences. First, the transition 
will be at a lower temperature citelee T c (ms) than in the stable state at T c 
by a factor -^/couplings. This follows by equating the temperature dependent 
mass tern in eq. [54] to zero at the critical temperature. 

Second, before the transition takes place, the system will undergo thermal 
inflation. Ej This happens in both states with complex and real values for the 
Wilson line in SO(10). 

Let us denote the VEV of the SO(10) Higgs H in the 16 representation by 
v. The relevant terms in the effective action are then: 

= l?H?\ H i H _ v 2f + T 2 (kn) 2 H^H + ... (35) 

Here k = 1, or ±~ according to wether p = — 1, or ±i. As long as the 
temperature is above the value T out = Tp- we have no VEV for the Higgs, 
and the energy density e of the vacuum equals ^(mv) 2 . This defines a second 
temperature scale T in = (mv)i . Below this scale the thermal contribution will 
become small with respect to the vacuum energy density till we reach the lower 
scale T out . Below that scale the symmetry breaks, the Higgs gets the VEV v 
and and the energy density becomes proportional to T 2 , from eq. |3^. So in 
between these temperatures we have inflation. 

There is a third scale, appearing in the Friedman equation coupling the 
radius R of a flat universe to its energy density: 

During the inflationary period the right hand side of this equation is 
approximately constant and equals the square of M v — . This scale is the 
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Hawking-Gibbons temperature Thg an d li es below T out . That is to say, it will 
not play any role in the inflationary period. This is a good feature, because 
Thg is a lower bound on the temperature during inflation. 



6 Conclusions 

We have given an overview of the use of the Wilson line in high temperature 
gauge field physics. 

The effective potential of the Wilson line does not only give rise to sym- 
metry breaking at high temperature. It provides in a natural way a gauge 
invariant infra red cut off and is useful in establishing the perturbative series 
for static quantities like free energy and Debije screening length. Many ap- 
plications are possible, e.g. Gallan-Symanzik type equations for the infrared 
behaviour of hot gauge theorycj. Very striking are its universal features: Long 
lived metastable states with potentialy very interesting consequences. 
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